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1.  Introduction 


In  support  of  the  U.S.  Anny  vision  for  increased  mobility,  survivability,  and  lethality,  the  U.S. 
Army  Research  Laboratory  (ARL)  has  recently  developed  the  forward-looking  ultra-wideband 
(UWB)  synthetic  aperture  radar  (SAR)  ( 1 ,  2,  3).  The  radar  transmits  and  receives  short  impulses 
that  cover  the  frequency  spectrum  from  300  MHz  to  3000  MHz.  Conventionally,  very  high 
speed  analog-to-digital  converters  (ADCs)  are  required  to  directly  digitize  these  time  domain 
impulses.  Although  this  conventional  technique  allows  the  radar  design  to  be  simple,  it  is  too 
costly  to  build  the  high  speed  ADCs  that  can  directly  digitize  the  wide-bandwidth  returned 
impulses.  Therefore,  ARL  has  developed  a  sampling  technique  that  allows  us  to  employ 
inexpensive  ADCs  to  digitize  wide  bandwidth  signals.  This  technique  is  called  Synchronous 
Impulse  REconstruction  (SIRE).  This  is  a  modified  and  enhanced  version  of  the  equivalent 
time-sampling  technique  that  is  widely  used  in  commercial  digital  storage  oscilloscopes  ( 4 )  and 
some  radar  (5).  As  mentioned,  the  advantage  of  this  technique  (or  other  equivalent  time¬ 
sampling  techniques)  is  that  it  relieves  the  speed  requirements  for  the  ADCs.  However,  the  data 
acquisition  time  for  signal  reconstruction  increases.  If  the  radar  is  stationary  during  this  data 
acquisition  cycle,  the  radar  signal  will  be  perfectly  reconstructed  from  the  returned  impulses. 
However,  this  assumption  is  not  true  since  the  relative  position  between  the  radar  and  the  targets 
is  no  longer  negligible  during  the  data  acquisition  cycle  because  of  the  motion  of  the  radar 
platform.  This  is  especially  true  for  forward-looking  radar  since  the  radar  moves  toward  the 
imaging  area.  This  relative  motion  between  the  radar  and  the  targets  during  the  data  acquisition 
cycle  results  in  severe  phase  and  shape  distortions  in  the  reconstructed  signal.  This  in  turn 
results  in  poor  focus  quality  and  low  signal-to-noise  level  in  SAR  imagery. 

This  report  describes  a  signal-processing  technique  that  removes  the  phase  and  shape  distortions 
from  the  radar  signal,  which  are  attributable  to  the  motion  of  the  platfonn.  This  technique  results 
in  SAR  imagery  with  significant  improvement  in  focus  quality  and  signal-to-noise  level.  This 
has  been  tested  and  verified  with  simulated  and  measured  radar  data.  This  technique  could  be 
applied  for  any  time-based  impulse  radar  system  that  experiences  the  relative  motion  between 
the  radar  and  the  targets  during  the  data  acquisition  cycle. 


2.  Signal-Processing  Technique  to  Compensate  for  Forward  Motion 


2.1  The  ARL  SIRE  Technique 

The  ARL  SIRE  radar  system  (3)  employs  an  Analog  Devices  12-bit  80-MHz  ADC  to  digitized 
returned  radar  signals.  However,  the  ADC  is  clocked  at  the  system  clock  rate  of  40  MHz.  From 
the  basic  sampling  theory  ( 4 ),  it  is  not  possible  to  reconstruct  the  wide-bandwidth  signal  (300 
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MHz  to  3000  MHz)  since  this  ADC  rate  is  much  slower  than  the  required  minimum  Nyquist1 
sampling  rate  (6000  MHz).  However,  by  using  the  synchronous  time  equivalence  sampling 
technique,  we  can  achieve  a  much  higher  equivalent  sampling  rate.  Figure  1  provides  a 
graphical  representation  of  the  SIRE  data  acquisition  technique. 


Figure  1.  The  ARL  synchronous  impulse  reconstruction  data  acquisition  scheme.  (This  is  a  modified  and 
enhanced  version  of  the  equivalent  time-sampling  technique.) 

The  ADC  sampling  period  is  A,;  the  value  of  this  parameter  in  figure  1  is  25  nsec,  which 
corresponds  to  an  analog-to-digital  (A/D)  sampling  rate  of  40  MHz.  The  number  of  samples  for 
each  range  profile  is  denoted  by  N,  which  is  equal  to  7  in  our  example.  This  corresponds  to  a 
range  swath  of  30  m.  The  system  pulse  repetition  frequency  ( PRF)  is  1  MHz.  The  system  pulse 

repetition  interval  ( PRI  =  — - —  ),  i.e.,  the  inverse  of  PRF,  is  1  micro-second  ((mu)s).  Each 

PRF 

aliased  (slowly  sampled  at  A/D  speed)  radar  record  is  measured  M  times  (1024  in  this  example) 
and  the  records  are  integrated  to  achieve  a  higher  signal-to-noise  level.  After  M  repeated 
measurements  of  the  same  range  profile  are  summed,  the  first  range  (fast-time)  bin  is  increased 
by  A.  Thus,  the  next  group  of  M  range  profiles  is  digitized  with  a  timing  offset  of  A  with  respect 
to  the  transmit  pulse.  The  parameter  A  represents  a  time  sample  spacing  that  satisfies  the 
Nyquist  criterion  for  the  transmitted  radar  signal.  For  the  cited  case  in  figure  1,  this  time  sample 
spacing  is  129.53  pico-seconds,  which  corresponds  to  a  sampling  rate  of  7.72  GHz.  This 


^Nyquist  is  the  person  who  invented  this  theory. 
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effective  sampling  rate  is  sufficient  for  the  wide-bandwidth  radar  signal  (300  MHz  to  3000 
MHz).  Note  that  the  number  of  interleaved  samples  is 


which  is  193  in  figure  1.  After  K  groups  of  M  pulses  are  transmitted  and  the  return  signals  are 
digitized  and  summed  by  the  Xilinx  Spartan2  3  field  programmable  gate-array  (FPGA),  this 
results  in  a  radar  record  of  N.K  samples  with  an  equivalent  of  fast  sample  spacing  of  A.  The  total 
time  to  complete  one  data  acquisition  cycle  is  N.K.PRI,  which  is  197.6  msec  in  this  case.  Please 
note  that  during  the  entire  data  acquisition  cycle  period  (197.6  msec),  the  relative  position 
between  the  radar  and  the  targets  is  assumed  to  be  stationary.  Table  1  summarizes  the 
parameters  used  by  the  SIRE  data  acquisition  technique. 

Table  1 .  Summary  of  radar  parameters. 


Radar  pulse  repetition  frequency  ( PRF) 

1  MHz 

Radar  pulse  repetition  interval  (PRI) 

le-6  sec 

ADC  sampling  rate 

40  MHz 

ADC  sampling  period 

25  nsec 

Number  of  ADC  (slow)  range  gates  (A) 

7 

Interleaving  factor  ( K) 

193 

Number  of  repeated  measurements  for  averaging  (M) 

1024 

Total  number  of  range  gates  (N.K) 

1351 

Effective  sampling  period  (time-equivalent) 

129.53e-12  sec 

Effective  sampling  rate  (time-equivalent) 

7.72  GHz 

Total  data  acquisition  time 

197.6  msec 

2.2  Perfect  Reconstruction  of  the  Radar  Signal  With  Stationary  Radar  Platform 

Let  us  consider  the  first  simulation  case.  The  radar  is  situated  10  m  away  from  a  point  target. 

The  radar  transmits  impulse  signals  to  the  point  target.  The  receiver  performs  the  data 
acquisition  on  returned  signals  using  the  sampling  technique  described  in  previous  section.  After 
M.K  (1024x193)  pulses  are  transmitted  and  received,  the  data  acquisition  cycle  is  completed  and 
the  radar  waveform  is  reconstructed.  In  this  case,  since  the  radar  is  stationary,  the  reconstructed 
waveform  is  perfect.  Figure  2  shows  both  the  time  domain  and  frequency  domain  of  the 
reconstructed  waveform  with  stationary  radar. 


2Xilinx  Spartan  is  a  registered  trademark  of  Xilinx,  Inc. 
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Time  Domain  of  Reconstructed  Frequency  Domain  of  Reconstructed 


Figure  2.  Perfect  reconstruction  with  stationary  radar. 

2.3  Phase  and  Amplitude  Distortions  in  Reconstructed  Radar  Signals  With  Moving 
Platform 

Next,  we  want  to  illustrate  the  problem  when  the  radar  is  moving  during  the  data  acquisition 
cycle.  This  is  the  same  simulation  case  as  in  previous  section  except  that  the  radar  is  moving 
toward  the  target  at  the  speed  of  5  miles  per  hour.  In  figure  3,  we  show  the  reconstructed 
waveform  of  the  moving  case  versus  the  stationary  case.  In  the  time  domain  plot  of  figure  3,  we 
can  notice  a  significant  phase  shift  in  the  reconstructed  waveform  compared  with  the  ideal  case 
(stationary).  The  phase  information  is  crucial  for  the  SAR  image  formation  process  since  radar 
signals  are  coherently  processed  by  the  image  former.  In  addition  to  the  phase  shift,  the 
frequency  domain  plot  of  the  moving  case  (figure  3)  indicates  that  there  is  a  distortion  in  the 
shape  of  the  reconstructed  waveform.  The  amount  of  phase  shift  and  shape  distortions  changes 
as  the  distance  from  the  radar  to  the  target  varies.  Figure  4  shows  the  reconstruction  wavefonn 
of  another  simulation  case  that  is  similar  to  the  case  in  figure  3,  except  that  the  distance  between 
the  radar  and  the  target  is  15  m  instead  of  10  m.  In  addition  to  the  distance  between  the  radar 
and  the  target,  the  speed  of  the  platform  contributes  to  the  amount  of  distortions.  Figure  5  shows 
the  reconstructed  wavefonn  similar  to  the  simulation  case  in  figure  4,  except  that  the  radar 
platform’s  speed  is  increased  to  12  miles  per  hour. 
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Figure  3.  Phase  and  shape  distortions  in  reconstructed  waveform  (red).  (The  distance  between  the  radar  and  the 
target  is  10  m.  The  speed  of  the  radar  platform  is  5  miles  per  hour.  The  ideal  reconstructed  waveform 
[stationary  radar]  is  shown  in  blue.) 
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Figure  4.  Phase  and  shape  distortions  in  reconstructed  waveform  (red).  (The  distance  between  the  radar  and  the 
target  is  15  m.  The  speed  of  the  radar  platform  is  5  miles  per  hour.  The  ideal  reconstructed 
waveform  [stationary  radar]  is  shown  in  blue.) 
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Figure  5.  Phase  and  shape  distortions  in  reconstructed  waveform  (red).  (The  distance  between  the  radar  and 
the  target  is  15  m.  The  speed  of  the  radar  platform  is  12  miles  per  hour.  The  ideal  reconstructed 
waveform  [stationary  radar]  is  shown  in  blue.) 


2.4  Technique  to  Remove  the  Phase  and  Shape  Distortions  From  Reconstructed 
Waveforms 

In  this  section,  we  discuss  the  technique  to  remove  the  phase  and  shape  distortions  from  the 
reconstructed  waveforms.  We  show  the  results  using  simulated  data  and  real  radar  data. 

The  phase  and  shape  distortions  in  the  reconstructed  waveform  can  be  explained  as  follow. 

Figure  6  shows  the  data  acquisition  scheme  when  the  radar  is  moving  during  the  data  acquisition 
cycle.  For  the  reconstructed  waveform,  the  effective  sampling  period  is  Am,  which  no  longer  has 
the  same  value  as  A  in  figure  1 .  After  transmitting  and  receiving  the  first  M  pulses,  the  radar 
starts  the  second  group  of  M  pulses  with  the  timing  offset  of  A  from  the  previous  group.  During 
the  data  acquisition  time  for  the  first  M  pulses,  in  this  case,  the  radar  has  traveled  a  distance  of 

d  =  V.M.PRI  (2) 

in  which  v  is  the  speed  of  the  radar  during  this  acquisition  group  of  M  pulses,  and  PR1  is  the 
pulse  repetition  rate  mentioned  in  section  2.1. 
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Figure  6.  The  ARL  synchronous  impulse  reconstruction  data  acquisition  scheme  when  the  radar  is  moving  during 
the  data  acquisition  cycle. 

The  distance  d  corresponds  to  the  elapsed  time  of 

2d_ 

c 

=  2vMPRl  (3) 

c 

o 

in  which  c  =  3  x  10  is  the  speed  of  light. 

Thus,  the  effective  sampling  period  for  the  reconstructed  waveform  is 

A„,  =  A  +  t 

m 

,  2  V.M.PRI 

=  A+ -  (4) 

c 

From  equation  4,  the  effective  sampling  period  for  the  reconstructed  waveform  is  varied  with  the 
radar’s  instantaneous  speed.  This  generates  the  phase  and  shape  distortion  in  the  reconstructed 
waveform.  Another  source  of  artifact  in  the  reconstructed  waveform  is  depicted  in  figure  6.  The 
last  samples  of  the  first  range  (measured  at  the  slower  A/D  rate)  overlap  with  the  beginning 
samples  of  the  second  range  gate.  This  creates  the  spurious  glitches  in  the  reconstructed 
waveforms  (shown  in  figures  4  and  5).  From  equations  2  to  4,  we  assume  that  the  radar  speed  v 
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is  constant  during  the  entire  data  acquisition  cycle.  This  is  not  a  bad  assumption  since  the  radar 
is  moving  slowly  and  its  speed  should  not  change  too  much  during  the  period  of  197  msec. 

Let  sm  be  the  reconstructed  waveform  with  phase  and  shape  distortions  attributable  to  the  radar 
platform’s  motion  during  the  data  acquisition  cycle.  Given  the  average  speed  of  the  radar 
platform  at  the  time  the  measurement  is  made,  we  want  to  compute  the  ideal  reconstruction  s 
from  the  distorted  sm  as  if  the  radar  is  stationary  during  the  data  acquisition  cycle.  The  system 
employs  a  differential  global  positioning  system  (GPS)  system  to  measure  the  radar  locations 
along  the  path.  In  addition,  the  GPS  time  stamps  are  recorded  with  the  radar  data  stream.  With 
the  radar  coordinates  and  time  stamps  infonnation  from  two  successive  locations,  the  average 
speed  of  the  radar  at  every  location  can  be  computed.  The  procedure  to  compensate  for  the 
phase  and  shape  distortions  in  the  radar  signatures  is  as  follows. 

First,  the  distorted  radar  signal  sm  is  measured. 

Next,  the  avarage  speed  at  this  position  is  computed  as 


in  which  (jc/_i,  y,-i,  z;-i)  are  the  coordinates  of  the  radar  at  the  previous  acquisition  cycle  (f-i), 
and  (xj,  Vi,  zt)  are  the  coordinates  of  the  radar  at  the  current  acquisition  cycle  (/,). 

The  equivalent  sampling  period  for  this  acquisition  cycle  Am  is  then  computed  as  specified  by 
equation  4.  The  distorted  radar  signal  sm(n)  is  divided  into  N  sub-signals  as 

smi  =  sm  (/ ) ,  where  1  +  (/  - 1) K  <  l  <iK  and  1  <  i  <  N .  (6) 

Each  sub-signal  smi  is  then  re-sampled  with  the  use  of  the  computed  sampling  period  Am 

(7) 

in  which  P  is  the  interpolation  process  that  re-samples  the  signal  from  the  sampling  period  Am  to 
the  new  sampling  period  Am. 

The  reconstructed  waveform  is  the  concatenation  of  N  sub-signals 

■s  =  k,-s1,-,-s„-,-sJV]  (8) 

2.5  Results 

2.5.1  Results  Using  Simulated  Data 

We  applied  the  technique  described  in  section  2.4  to  the  distorted  signature  (red  in  figure  4).  The 
time  and  frequency  results  are  plotted  in  figure  7.  The  distorted  signature  is  perfectly 
reconstructed.  The  shape  distortion  is  removed  and  its  phase  is  aligned  with  the  original 
stationary  signature. 
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Figure  7.  The  technique  is  applied  to  the  simulated  data  (in  figure  4).  (This  figure  shows  the  perfect 
reconstruction  with  the  processing  technique.) 

2.5.2  Results  Using  Measured  Radar  Data 

The  ARL  SIRE  radar  system  has  an  array  of  16  receivers  that  form  a  2-m  aperture  and  two 
transmitters  that  are  situated  at  the  ends  and  above  the  receiving  array.  The  radar  transmits  M 
impulses  using  the  left  transmitter  and  simultaneously  digitizes  and  integrates  the  return  radar 
signals  from  16  receiving  channels.  The  radar  then  switches  to  the  right  transmitter  and  repeats 
the  same  process.  The  integration  is  perfonned  by  the  FPGA  in  the  radar  hardware.  Each  time  a 
transmitter  is  selected,  the  radar  generates  one  data  frame  that  includes  data  from  all  16 
receivers.  Figure  8  shows  seven  frames  of  radar  signatures  from  an  8-foot  trihedral  as  the  radar 
moves  toward  the  target  at  the  speed  of  5  mph.  In  each  data  frame,  we  should  expect  two 
responses:  a)  the  first  large  response  from  the  trihedral,  and  b)  a  much  weaker  response  because 
of  the  transmitter  artifact.  However,  we  can  notice  the  obvious  shape  distortion  that  occurred  in 
frames  1,  2,  4,  and  6.  In  frames  2  and  6,  the  main  responses  split  while  in  frames  1  and  4,  the 
secondary  responses  split.  In  addition  to  the  shape  distortion,  the  phase  errors  are  present.  We 
applied  the  technique  described  in  section  2.4  to  these  data.  Figure  9  shows  the  data  after 
processing.  Notice  that  the  phase  and  shape  distortions  of  the  trihedral  signatures  are  corrected 
in  figure  9.  Figure  10  shows  the  time  domain  response  of  the  same  trihedral  before  and  after 
forward  motion  processing. 


9 


Figure  8.  Raw  radar  data  measured  by  ARL  SIRE.  (The  radar  moves  at  5  mph  toward  an  8-ft  trihedral. 

This  figure  shows  the  distortions  in  radar  signatures  when  the  radar  is  moving  during  the  data 
acquisition  cycle.) 


Figure  9.  Data  in  figure  8  have  been  processed  with  the  described  processing  technique.  (The  distortions  in  radar 
signatures  have  been  removed.) 
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Figure  10.  A  down-range  profile  through  a  trihedral  before  and  after  processing. 


We  have  seen  the  effect  of  radar  forward  motion  on  the  radar  signature.  After  all  the  signal¬ 
processing  algorithms  are  perfonned  on  radar  phase  history  data,  the  back-projection  image 
formation  algorithm  is  applied  to  the  radar  data  to  produce  a  SAR  image. 

Over  the  last  decade,  ARL  has  been  using  the  back-projection  image  formation  for  various 
applications  (8  through  11).  In  this  forward  imaging  case,  the  back-projection  algorithm  uses  a 
two-dimensional  (2-D)  aperture  to  perform  the  imaging  process  that  coherently  integrates  data 
from  this  2-D  aperture.  In  this  case,  without  compensation,  the  phase  and  shape  distortions  in 
the  original  radar  signatures  would  result  in  poorly  focused  SAR  imagery.  Figure  1  la  (left) 
shows  a  SAR  image  with  data  from  the  forward-looking  ARL  SIRE  radar.  The  SAR  image 
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covers  an  area  of  25  m  in  cross-range  direction  and  64  m  in  down-range  (or  forward)  direction. 
In  this  run,  the  radar  moves  very  slow  at  an  average  speed  of  1  m/sec.  All  the  signal  processing 
steps  are  applied  to  the  data  except  this  forward  motion  compensation  processing.  Figure  lib 
(right)  shows  the  same  SAR  image  (plotted  with  the  same  dynamic  range)  with  all  the  signal¬ 
processing  steps  (that  include  the  forward  motion  compensation  processing)  applied  to  the  radar 
data.  It  is  obvious  that  without  the  forward  motion  processing,  the  targets  in  figure  1  la  do  not 
preserve  the  original  shapes  and  significantly  suffer  the  loss  of  their  radar  cross-sectional  values. 
Also  note  that  the  details  of  the  clutter  strip  on  the  right  side  of  the  SAR  image  are  lost  in  figure 
11a.  Figure  12  shows  the  image  chips  of  four  in-scene  targets  with  (right  column)  and  without 
(left  column)  the  forward  motion  processing.  The  forward  motion  processing  improves  the 
peak’s  levels  of  the  targets  from  5  dB  to  13  dB  in  this  case.  More  importantly,  without 
processing,  the  shape  distortions  will  lead  to  the  wrong  frequency  contents.  For  a  typical  target 
detection  system  (11),  the  textures  and  frequency  contents  of  the  targets  play  important  roles  in 
the  discrimination  of  targets  against  other  confuser  classes. 


Without  Forward  Motion  Processing  With  F orward  Motion  Processing 


-10  -5  0  5  10  -10  -5  0  5  10 


_ Cross-Range  (meter) _ Cross-Range  (meter) 

Figure  11.  SAR  image  using  data  from  the  forward-looking  ARL  SIRE  radar.  (The  SAR  image  on  the 

left  does  not  have  the  application  of  the  forward  motion-processing  algorithm.  The  SAR  image 
on  the  right  shows  the  result  after  the  forward  motion-processing  technique  has  been  applied 
to  the  radar  data.) 
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Without  Forward  Motion 
Processing 


110.6  dB 


With  F  orward  Motion 
Processing 


90.4  dB 


Figure  12.  SAR  image  chips  of  targets  before  (left)  and  after  (right)  the  processing 

of  the  forward  motion  algorithm.  (The  forward  motion  processing  improves 
the  peak’s  levels  of  the  targets  from  5  dB  to  13  dB  in  this  case.) 
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3.  Summary 


ARL  has  designed  and  developed  the  SIRE  UWB  SAR  radar  in  support  of  the  U.S.  Anny  vision 
for  increased  mobility,  survivability,  and  lethality.  The  radar  is  based  on  time-domain  wideband 
impulses.  For  this  radar,  ARL  designed  and  implemented  a  data  acquisition  technique  called 
SIRE  that  allowed  us  to  employ  relatively  slow  ADC  (40  MHz)  to  digitize  wideband  signals 
(>3000  MHz).  However,  the  scheme  assumed  that  the  radar  and  targets  are  stationary  during  the 
data  acquisition  cycle  when  in  reality,  the  target  signatures  did  suffer  the  distortions  in  phase  and 
shape  because  of  the  radar  motion.  The  phase  error  would  lead  to  significant  loss  in  target  radar 
cross-sectional  values  in  resulting  SAR  imagery.  The  shape  errors  would  destroy  the  frequency 
contents  of  the  targets  and  thus  the  ability  to  discriminate  targets  from  other  confuser  classes. 

This  report  described  a  signal-processing  method  that  we  designed  to  recover  the  accuracy  of  the 
target  signatures  that  were  affected  by  the  radar  motion.  We  applied  the  signal-processing 
method  to  simulated  data  and  measured  data  from  ARL  SIRE  radar.  With  the  simulated  data,  we 
showed  the  perfect  reconstruction  using  this  method.  With  the  measured  radar  data,  we 
demonstrated  significant  improvement  in  the  resulting  SAR  image.  The  radar  cross  section  of 
targets  improved  from  5  dB  to  13  dB.  The  correct  shapes  of  target  signatures  are  preserved. 
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